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Abstract 

Simple exclusion processes for particles moving along two parallel lattices and jumping 
between them are theoretically investigated for asymmetric rates of transition between the 
channels. An approximate theoretical approach, that describes the particle dynamics ex- 
actly in any vertical cluster of two parallel sites and neglects the correlations between the 
different vertical clusters, is applied to calculate stationary-state density profiles, currents 
and phase diagrams. Surprisingly, it is found that asymmetry in the coupling between the 
channels leads to a very complex phase behavior that is very different from two-channel sim- 
ple exclusion processes with symmetric coupling. There are seven stationary-state phases 
in the simple exclusion processes with asymmetric transition rates between the channels, 
in contrast to three phases found for the systems with symmetric coupling. In addition, 
a new maximal-current phase with a domain wall in the middle of the lattices, that has 
no analogs in other exclusion processes, is observed. Although the explicit calculations are 
presented only for the case of full asymmetry, when the particles can only jump between 
the channels in one direction, the properties of two-channel simple exclusion systems with 
general asymmetry are also discussed. Theoretical predictions are in excellent agreement 
with extensive computer Monte Carlo simulations. 
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1 Introduction 



Asymmetric simple exclusion processes (ASEP) play a critical role for understanding multiple 
non-equilibrium phenomena in chemistry, physics and biology mi2]- ASEPs have been extensively 
studied and applied for description of kinetics of biopolymerization jSj, protein synthesis [HE], 
transport of motor proteins in biological cells [Hj, polymer dynamics in dense medium car 
traffic processes [8 and modeling of ant trails 

The majority of studied asymmetric exclusion models analyze the multi-particle dynamics 
along a single lane d 12] • At the same time, the more realistic description of many processes, 
such as car traffic and the biological transport of motor proteins 13 E], suggest that it is important 
to study multi-lane ASEPs. There have been several theoretical investigations of two-lane simple 
asymmetric exclusion processes [TIH ITTl IT^ ITBl I14j . We introduced earlier two-channel simple 
exclusion models where particles can move along the channels and between them ^T]. Using a 
vertical cluster mean-field approach the stationary properties have been analyzed and compared 
with extensive computer Monte Carlo simulations. Only the case of symmetric coupling, when the 
particles jump with equal rates between the lanes, has been considered. Mitsudo and Hayakawa 
[T!^ extended these models to more general asymmetric rates of moving between the channels, and 
it was shown that kinks, or domain walls, in both channels synchronize their motion. However, 
the theoretical analysis has been based on the mean-field decoupling approximation that totally 
neglects all correlations in the system, and significant deviations from Monte Carlo computer 
simulations results have been observed, especially when the coupling is asymmetric. The effect 
of disorder in the transition rates and current-density relations for limited range of parameters 
have also been studied for two- lane exclusion processes with asymmetric coupling ,14. 

The goal of this paper is to investigate the general stationary-state properties of two-channel 
exclusion processes with asymmetric coupling for all possible sets of parameters. We consider 
the effect of asymmetry in the transition rates between the lanes by using the vertical cluster 
mean-field approach, that takes into account the correlations between the channels. In addi- 
tion, extensive computer Monte Carlo simulations are performed in order to test the theoretical 
predictions. 

The paper is organized as follows. Theoretical description of the model and mean-field analy- 
sis are presented in Section 2. The results of Monte Carlo computer simulations and comparison 
with theoretical predictions are discussed in Section 3. Summary and conclusions are presented 
in the final Section 4. 

2 Theoretical Description 
2.1 Model 

We consider a system of two parallel one-dimensional lattices where identical particles can move 
along the channels and between them, as shown in Fig. 1. Each lattice has L sites, and every 
site can be occupied by no more than one particle or it can be empty. At every time step a 
site is randomly chosen from the lane 1 or 2. In the bulk of the system the particle dynamics 
is described by the following rules. A particle can change the lane with the rate wi or W2 from 
the channel 1 or 2, correspondingly, if the vertical neighboring site is available - see Fig. 1. The 
particle at site 1 < i < L can also move from left to right along the same channel to site i + 1, if 
this site is empty. The rates for horizontal moves from the site i depend on the occupancy state 
of the site i on another lane. The particle on lane 1 moves to the right with the rate 1 — wi if 
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the vertical neighboring site is not occupied, otherwise it jumps with the rate 1. Similarly, for 
the particle on lane 2 the horizontal transition rates are equal to 1 — W2 (the vertical neighbor is 
empty) or 1. These rules satisfy the condition that the total probability per unit time of leaving 
the site i (in any direction) is always equal to 1 [TTj . 

In addition, there are special entrance and exit dynamic rules at the boundaries. Particles 
can enter the system with the rate a if any of the first sites at each lane is not occupied. When a 
particle reaches the exit site L it can leave the system with the rate j3 if the exit site on another 
lane is occupied. Otherwise, the exit rates are /?(1 — Wi) and /5(1 — W2) for the channel 1 and 2, 
corresp ondingly. 
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Figure 1: Schematic view of the model for a two-channel ASEP with asymmetric coupling. 
Allowed transitions are shown by arrows. The inter-channel transition rates are equal to Wi and 
W2. The horizontal transition rates are 1 — Wi and I — W2 for the channel 1 and 2, correspondingly, 
if the vertical neighboring sites are empty. Otherwise, particle move with the rate 1. Entrance 
rates are equal to a at both lanes. Exit rates are equal to /5(1 — wi) and /3(1 — W2) for the lane 
1 and 2, respectively, if the exit vertical cluster is half-occupied. For the fully-occupied vertical 
cluster at the exit the rate of leaving is f3. 



For simplicity, in the present work we consider only the case of full asymmetry in the vertical 
transition rates with wi = 1 and W2 = 0. However, our analysis can also be straightforwardly 
applied for more general asymmetric couplings in two-channel ASEPs. Note also that there is a 
particle-hole symmetry in the system, and it can be mapped into itself by exchanging everywhere 
the labels 1 and 2 and the rates a and (3. This symmetry is important for understanding the 
stationary-state properties of two-lanes exclusion processes. 



2.2 Mean- field Analysis 

Theoretical study of two-channel ASEP with symmetric coupling JT] indicates that correlations 
between the channels strongly influence the stationary-state properties of the system, however the 
horizontal correlations inside the lanes are relatively weak. The effect of inter-channel correlations 
is even larger for two-channel ASEPs with asymmetric coupling [Ij . This suggests that a cluster 
mean-field approach jTT], that explicitly takes into account the correlations inside the vertical 
cluster of lattice sites, is the most appropriate and convenient theoretical tool for analyzing these 
systems. 

The basic quantities of the vertical cluster mean-field approach are the probabilities to find 
any vertical cluster in one of four possible states, as shown in Fig. 2. Considering the lattice 
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sites far away from the boundaries of the system, it is assumed that the occupation of vertical 
clusters is independent of the position along the channels. We define Pn as a probability to find 
a vertical cluster with both lattice sites filled, Pio and Pqi as probabilities to have a half-empty 
vertical cluster with a particle at the channel 1 or 2, respectively, and Pqo as a probability to 
have no particles at both lattice sites. The conservation of probability requires that 

Pn + Poi + Pio + Poo = 1. (1) 

In addition, the bulk densities at each channel can be expressed in terms of the vertical cluster 
probabilities, 

pi = Pn + Pio, P2 = Pn + Poi- (2) 
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Figure 2: Four different configurations for vertical clusters of lattice sites. Pn, Pio, Pqi and Pqo 
are the corresponding probabilities for each state. 

The dynamics of the system can be described via Master equations for evolution of every 
vertical cluster state. Specifically, for the fully occupied cluster we have 



dp 



11 



dt 



PiiPio + PiiPoi + PqiPio ~ 2PiiPoo — PnPio — PnPoi. 



(3) 



In the limit of t ^ cxo the system reaches a stationary state with = 0, and this equation 
simplifies into 

PoiPio = 2PnPoo. (4) 
Similarly, for the half- filled vertical clusters it can be shown that 
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PllPoo + PnP 
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that at the stationary-state limit = 0) reduces to 

-Plo(1 ~ -Pii) = 2PiiPoo- 



(5) 



(6) 



From Eqs. (HJ, (jH) and © it is possible to calculate the stationary-state probabilities of the differ- 
ent vertical cluster states, and consequently all properties of two-channel ASEPs with asymmetric 
coupling can be obtained. 

In the bulk of the system the stationary currents are given by the following expressions. 



J, 



bulk,l 



-Pii(l — -Pii — Pio), Jbuik,2 — {Pn + -Poi)(l ~ -Pii ~ -Poi)- 



(7) 



The particle currents at the boundaries are different form the bulk expressions. For entrance we 
obtain 

Jentr,! = ~ Pn — Pio) , Jentr,2 = «(1 — Pn — Pm)] (§) 
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while at the exit it can be shown that 



Jexit,l — PPn, Jexit,2 — P{Pll + Pqi) ■ (9) 

Ahhough the currents on the individual lanes can differ from each other, the overall current of 
the system at the stationary state is always constant, 

Jtotal = Jbulk,l + Jbulk,2 = Jentr,l + Jentr,2 = Jexit,! + Jexit,2- (10) 

Solving together Eqs. (PJ), Q and © leads to a conclusion that at large times there are three 
possible cases. In one region of phase space for all bulk lattice sites we have 

Pio = Pii = 0, (11) 

which means that the bulk density at the lane 1 is zero. In another region of the parameter's 
space 

Pio = Poo = 0, (12) 

that leads to p2 = 1, i-e., the channel 2 is fully occupied by the particles. There is also a case 
when the two-channel system can be divided in two parts, with Eq. (fTT|) valid in one part while 
Eq. (jl2j) is satisfied in another part. 

Let us consider first the region in the parameter's space when the system is described by Eq. 
(|TT|). In this case there are only two possible states for the vertical clusters, the cluster with both 
empty sites (00) and the cluster with the occupied site on lane 2 (01). The expressions for the 
particle currents are significantly simplified, 

Jbulk,l = 0, Jbulk,2 = Poi(l ~ Poi), 
Jentr,l = Oi, Jentr,2 = Q^(l — Poi) 

Jexit,l = 0, Jexit,2 = PPol- (13) 

Then the dynamics of the two-lane system can be viewed as an effective one-channel transport 
with "particles" given by (01) vertical clusters and "holes" being (00) vertical clusters. The 
particles are entering into the system with an effective rate aeff 7^ a, that can be calculated 
from Eqs. (jl3p . This is due to the fact that there are two entrance currents. At the same time 
the effective particles are exiting with the rate (3 since there is only one exit current. It is known 
m 12] that the corresponding one-channel ASEP has low-density, high-density and maximal- 
current phases, and we conclude that there are also three phases, labeled as (0,LD), (0,HD) and 
(0,MC), in this part of the phase space for the two-channel ASEP with asymmetric coupling. 
The indexes in these labels for the different phases reflect the fact that pi = is always in the 
bulk of the system, while the state of the channel 2 depends on the entrance and exit rates. 

The phase boundaries and the effective entrance rate aeff can be determined in the following 
way. In the low-density phase (0,LD) the entrance current determines the total current in the 
system, namely, 

Jentr,! + Jentr,2 = Jbulk,! + Jbulk,2- (14) 

From Eqs. (fTSj) we have 

« + a(l - Poi) = ae//(l - Poi) = Poi(l - ^i), (15) 
with Pqi = aeff. It can be shown that 



1 + a- J{l + ay -8a 1 - a + J{1 + a)^ - 8a 
aeff = Poi = p2 = , Poo = • (16) 
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This phase exists for Oe// < /? and Og// < 1/2, that leads to the following conditions, 



P> ^ , a<l/6, (17) 



with the total current in this phase given by 

Jtotal = ae//(l - «e//) = (tt/2) 



3-a + J{l + ay -8a 



(18) 



In the phase (0,HD) the exit processes determine the particle dynamics in the channel 2, and 
this phase exists for f3 < aeff and f3 < 1/2. It can be easily calculated that in this phase 

Poi=P2 = l-/5, Poo = 13, Jtotal = 13(1-13). (19) 

The conditions aeff > 1/2 and (3 > 1/2 specify the phase (0,MD) with the maximal current at 
the channel 2. The stationary-state properties of this phase are given by 

POI = Poo = 1/2, Jtotal = 1/4. (20) 

It is also important to note that because the total particle current through the system cannot 
be larger than 1/4 all three phases that satisfy the Eq. (fTT|) cannot exist for a > 1/2. 

Similar calculations can be performed for the region of the phase space where Eq. ()12|1 is 
valid. Since in this case only the fully filled (11) and half-filed (01) vertical clusters can exist in 
the bulk, that leads to p2 = 1 and the system again can be mapped into effective one-channel 
ASEP. Then there are three possible phases, called (LD,1), (HD,1) and (MC,1), can be found 
for this range of parameters. The particle currents are given by 



Jbulk,l — -Pll(l 


-Pn), 


Jbulk,2 


= 0, 


JentrA = a{l 


-Pu), 


Jentr,2 


= 0, 


J exit, 1 ~ 


= PPu, 


J exit,2 


= p. 



(21) 

The (11) vertical cluster play a role of new effective particles that enter the system with the rate 
a and exit with an effective rate /3e//- From Eqs. ^IT\ it can be shown that 



Peff - ^ • \2Z) 

The phase (LD,1) exists for a < (3eff and a < 1/2. For this phase we obtain 

Pu=pi = a, Poi = l-", Jtotal = a{l - a). (23) 
The conditions for the (IID,1) phase are given by Peff < « and Peff < 1/2, that leads to 



1 + /3 - V(l + /3)2 - 8/3 
a> ^, (3<l/6. (24) 



In this phase the stationary-state current is equal to 



Jtotal = /3e//(l - f3eff) = {(3/2) [s - /3 + ^{1 + ^)^-8(3 



(25) 
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while the densities are 



Pii = Pi 



01 



(26) 



The third possible phase is (MC,1), defined for f3eff > 1/2 and a > 1/2. In this phase the 
steady-state properties are given by 



01 



Pii = 1/2, Ji 



total 



1/4. 



(27) 



Using the same arguments as above, all three phases described by Eq. (jl2j) cannot be found for 
/5 > 1/2. Note that these results can also be obtained by utilizing the particle-hole symmetry of 
the two-channel ASEPs with asymmetric coupling. The symmetry arguments also suggest that 
P is a phase boundary between (0,HD) and (LD,1) phases. 



a 



Finally, for a > 1/2 and /3 > 1/2 we have a situation when in one part of the two-channel 
system Eq. flllll is valid, while in the other part Eq. (fT^ determines the stationary behavior. 
We call this phase (MC,MC), and the total current in the system is equal to 1/4. Because of 
the particle-hole symmetry the boundary between two different parts is expected to be found 
exactly at the middle of the lanes. 

Thus, the vector cluster mean-field analysis suggests that there are seven stationary phases 
in two-channel exclusion processes with asymmetric coupling. The predicted phase diagram is 
shown in Fig. 3. There are two types of phase transitions can be observed in the system. 
The first-order phase transitions (shown by solid lines in Fig. 3) involve a jump in the particle 
densities, while in the continuum phase transitions (dashed lines in Fig. 3) there are smooth 
changes in the density profiles. 




Figure 3: Phase diagram for the two-channel ASEP with full asymmetric inter-channel coupling 
{wi = 1 and W2 = 0). Lines correspond to theoretical predictions, symbols are from the computer 
Monte Carlo simulations. Solid lines describe the first-order phase transitions, while dashed lines 
represent the continuum phase transitions. 

Although our calculations are performed only for the case of full asymmetry in inter-channel 
coupling, similar calculations can be made for systems with general asymmetry {wi > W2 > 0) . 
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It is expected to also find seven stationary phases in the phase diagram, although the phase 
boundaries, currents and the density profiles will depend explicitly on the vertical transition 
rates. 

3 Monte-Carlo Simulations and Discussions 

Our theoretical approach based on the vertical cluster mean-field description of two-channel 
ASEP with asymmetric coupling predicts a complex phase diagram with many non-equilibrium 
phase transitions. It also allows to calculate stationary currents and density profiles in both 
channels. In order to test our theoretical results, a series of extensive Monte Carlo computer 
simulations have been carried out. 

Since the predictions of the mean-field treatment are valid only in the thermodynamic limit, 
i.e., L oo, our simulations has been performed for different lattice sizes with L = 100, 500, 
1000 and 5000. Most of our computer simulations results are obtained for L = 1000 for which 
it was found that finite-size effects can be neglected. For faster computations the so-called BKL 
algorithm [15j has been used. This method skips uneventful Monte-Carlo steps and it makes a 
corresponding correction in the calculation of the effective Monte Carlo time. Employing the 
BKL algorithm has been especially useful for the cases of phase coexistence and for slow dynamics 
for weak entrance and exit rates. The density profiles and currents in our simulations have been 
computed by averaging 10^ Monte Carlo steps per site, although at the phase boundaries to 
obtain more accurate description we used 10^ — 10^ steps per site. In order to be confident that 
the system reached the steady state, typically first 3-5 % of all Monte Carlo steps have been 
neglected. 

The phase diagram computed from Monte Carlo simulations is presented in Fig. 3, and it can 
be seen that our theoretical predictions are in excellent agreement with computer simulations 
results. More information can be extracted from the density profiles that are shown in Fig. 4. 
Theoretically calculated density profiles in both channels also agree quite well with the results 
from computer simulations, although there are deviations near the boundaries of the channels. 
This is due to the fact that in our approximate theoretical treatment it was assumed that 
probabilities of vertical clusters are independent of the position along the lattice. Obviously, this 
assumption is valid only in the bulk of the system far away from the boundaries. It should be 
noted that because of the mapping of two-channel ASEP with asymmetric coupling into effective 
one-channel exclusion process the densities near the entrance and exit can be calculated 

The analysis of Figs. 4a, 4b and 4c allows to understand the nature of (0,LD), (0,IID) and 
(0,MC) stationary phases. In all three phases particles enter into the both channels (note that 
Pi 7^ at the left end of the system), but they exit only from the lane 2. Because of asymmetric 
inter-channel coupling the particles on the lane 2 have enough time to escape into the lane 1, and 
most of the particle transport in the system is taking place only along the channel 2. The density 
profiles for (LD,1), (IID,1) and (MC,1) phases can be easily obtained from (0,IID), (0,LD) and 
(0,MC) densities by using the particle-hole symmetry arguments. In these three phases particles 
enter only through the channel 1 but there are two exiting currents. Because of the slow exiting 
processes and asymmetric coupling at large times the channel 2 is fully occupied by the particles, 
and the most of the particles in the system move mostly along the channel 1. 

The most surprising result of our theoretical calculations is the existence of (MC,MC) phase 
shown in Fig. 4d. It has a domain wall, positioned in the middle of the system, that separates 
two different density profiles. In this phase particles enter the system only via the channel 1 
but exit only through the channel 2, and there is an inter-channel flux of the particles near the 
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Figure 4: Density profiles for a) (0,LD) phase with a = 0.05 and /? = 0.8; b) (0,HD) phase 
with a = 0.25 and /3 = 0.3; c) (0,MC) phase a = 0.4 and (3 = 0.9; and d) (MC,MC) phase 
with a — P — 0.9. Symbols correspond to Monte Carlo computer simulations, fines describe the 
density profiles in each channel. Error bars, determined from standard deviations for simulations, 
are smaller than the size of the symbols. 
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Figure 5: Density profiles for pliase transition line between (0,LD) and (HD,1) phase with 
a = (3 = 0.2. Symbols are from Monte Carlo computer simulations, lines are from theoreti- 
cal calculations. 




Figure 6: Stationary particle currents for (5 = 0.6. Symbols are from Monte Carlo computer 
simulations, lines are from theoretical calculations. 

middle position. The physical origin for this phase is the fact that for a > 1/2 and (3 > 1/2 
the phases (0,MC) and (MC,1) can coexist together. Because of the particle-hole symmetry the 
domain wall that separates two phases is exactly in the middle of the system. To the best of our 
knowledge, there are no such phases observed in any other driven diffusive systems jH |2] . 

Our theoretical calculations also agree well with computer simulations for the density profiles 
at phase transitions as shown in Fig. 5. Linear density profiles shown in Fig. 5 correspond 
to a first-order phase transition between (0,LD) and (HD,1) phases. The overall stationary 
dynamics and phase boundaries can also be understood by using the phenomenological domain 
wall description [T^ . 

The dynamics of two-channel ASEP with asymmetric coupling can be well understood by 
analyzing the stationary currents presented in Fig. 6. Again our theoretical predictions match 
quite well the results from Monte Carlo computer simulations. For the fixed value of = 0.6 
and for a < 1/Q the system is in (0,LD) phase and both entrance currents are increasing. At the 
phase transition (a = 1/6) between (0,LD) and (0,MC) phase Jentr,2 reaches the maximum and 
it starts to decrease for larger values of a, while Jentr,i is still growing. At a = 0.5 the second 
entrance current disappears and the first entrance current reaches a constant value of 1/4, and 
the system enters into (MC,MC) phase. 



10 



4 Summary and Conclusions 



Two-lane simple exclusion processes with asymmetric inter-channel coupling have been inves- 
tigated. To account for correlations between the channels, we utilized a mean-field method 
that allows to calculate exactly the probabilities of vertical clusters but neglects the correlations 
along the channels. It is shown that the asymmetry in inter-channel coupling produces a complex 
stationary-state behavior. There are seven stationary phases, and the particle dynamics in all 
phases can be understood in terms of effective one-channel ASEPs. The new (MC,MC) phase 
displays a domain wall between two density profiles in the middle of the system, and it has no 
analogs in other simple exclusion models. The specific position of phase boundaries depend on 
the degree of asymmetry, and for symmetric vertical transition rates the phase diagram simplifies 
with only three possible phases [TT] . 

There are several extensions of this model that will be important to investigate. In this 
paper the homogeneous vertical transitions rates have been considered. It will be interesting 
to understand the stationary properties of two-channel ASEP with inhomogeneous coupling. 
Another extension is to study the effect of adding the equilibrium Langmuir kinetics process of 
association and dissociation to one or to both lattices. Understanding these processes can help to 
describe better complex low-dimensional transport phenomena in chemistry, physics and biology. 
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